A LOWER BOUND FOR NODAL COUNT ON DISCRETE AND METRIC 

GRAPHS 



GREGORY BERKOLAIKO 



Abstract. We study the number of nodal domains (maximal connected regions on which a function 
has constant sign) of the eigenfunctions of Schrodinger operators on graphs. Under certain genericity 
condition, we show that the number of nodal domains of the n-th eigenfunction is bounded below by 
t*"** , n — £, where £ is the number of links that distinguish the graph from a tree. 

Our results apply to operators on both discrete (combinatorial) and metric (quantum) graphs. 
They complement already known analogues of a result by Courant who proved the upper bound n 
for the number of nodal domains. 

To illustrate that the genericity condition is essential we show that if it is dropped, the nodal count 

• can fall arbitrarily far below the number of the corresponding eigenfunction. 

\ In the appendix we review the proof of the case I = on metric trees which has been obtained by 

CO ■ other authors. 

<n : 

> ■ 

\q • 1. Introduction 

cn : 

According to a well-know theorem by Sturm, the zeros of the n-th eigenfunction of a vibrating 
^ string divide the string into n "nodal intervals". The Courant nodal line theorem carries over one 

• half of Sturm's theorem to the theory of membranes: Courant proved that the n-th eigenfunction 
O ■ cannot have more than n domains. He also provided an example showing that no non-trivial lower 
^ ! bound for the number nodal domains can be hoped for in M d , d > 2. 

But what can be said about the number of nodal domain on graphs? Earliest research on graphs 

• concentrated on Laplace and Schrodinger operators on discrete (combinatorial) graphs. The functions 
c3 ■ on discrete graphs take values on vertices of the graph and the Schrodinger operator is defined by 

(Hiff) v = - V" v« + Qui>u, 

> ■ 

^ \ where the sum is taken over all vertices adjacent to the vertex u. 

o3 ' Gantmacher and Krein [13] proved than on a chain graph (a tree with no branching which can be 
thought of as a discretization of the interval) an analogue of Sturm's result holds: the n-th eigenvector 
changes sign exactly n — 1 times. But for non-trivial graphs the situation departs dramatically from 
its R d analogue. First of all, Courant 's upper bound does not always hold. There is a correction 
due to multiplicity of the n-th eigenvalue and the upper bound becomes^ [8] n + m — 1, where m 
is the multiplicity. In this paper we discuss another striking difference. If the number of cycles of 
a graph is not large, the graph behaves "almost" like a string: for a typical eigenvector, there is a 
lower bound on the number of nodal domains. 

To be more precise, let I be the minimal number of edges of the graph that distinguish it from a 
tree (a graph with no loops). In terms of the number of vertices V and the number of edges E, the 
number i can be expressed as £ — E — V + l. We show that, for a typical eigenvector, the number 
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1 We are talking here about the so-called "strong nodal domains" — maximal connected components on which the 
eigenfunction has a constant well-defined (i.e. not zero) sign 
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of nodal domains is greater or equal to n — £. In particular, on trees (£ = 0) the nodal counting is 
exact: the n-th eigenfunction has exactly n domains. Here by a "typical" eigenvector we mean an 
eigenvector which corresponds to a simple eigenvalue^! and which is not zero on any of the vertices. 
This property is stable with respect to small perturbations of the potential {q u }- 

Another graph model on which the question of nodal domains is well-defined is the so-called 
quantum or metric graphs. These are graphs with edges parameterized by the distance to a pre- 
selected start vertex. The functions now live on the edges of the graph and are required to satisfy 
matching conditions on the vertices of the graph. The Laplacian in this case is the standard 1- 
dimensional Laplacian. A good review of the history of quantum graphs and some of their applications 
can be found in |19|. 



The ideas that the zeros of the eigenfunctions on the metric trees behave similarly to the 1- 
dimensional case have been around for some time. Al-Obeid, Pokornyi and Pryadiev [U [231 [22] 
showed that for a metric tree in a "general position" (which is roughly equivalent to our genericity 
assumption [21 see Section [3]) the number of the nodal domains of n-th eigenfunction is equal to 
n. This result was rediscovered by Schapotschnikow [24| who was motivated by the recent interest 
towards nodal domains in the physics community [H [151 E] • 

Our result on the lower bound extends to the quantum graphs as well. Similarly to the discrete 
case, we prove that even for graphs with £ > 0, n — £ is a lower bound on the number of nodal 
domains of the n-th eigenfunction. 

The article is structured as follows. In Section [2] we explain the models we are considering, 
formulate our result and review the previous results on the nodal counting on graphs. The case of 
the metric trees has been treated before in [221 I2"lj . In the three remaining cases, metric graphs 
with £ > 0, discrete trees and discrete graphs with £ > 0, we believe our results to be previously 
unknown and in Section [3] we provide complete proofs. For completeness, we also include a sketch 
of the general idea behind the proofs of [221 121] in the Appendix. Finally, in the last subsection of 
Section [3] we show that when a graph does not satisfy our genericity conditions, the nodal count can 
fall arbitrarily far below the number of the corresponding eigenfunction. 



2.1. Basic definitions from the graph theory. Let G be a finite graph. We will denote by V the 
set of its vertices and by £ the set of undirected edges of the graph. If there exists an edge connecting 
two vertices v\ and we say that the vertices are adjacent and denote it by V\ ~ v%. We will assume 
that G is connected. 

Definition 2.1. A graph G is connected if for any Vi,v 2 £ V there is a sequence of distinct vertices 
«i, . . . u n leading from Vi to v 2 («i = v i, u n = v 2 and Uj ~ for j — 1, . . . n — 1). A graph G is a 
tree if for any v i and v 2 the sequence of Uj connecting them is unique. 

The number of edges emanating from a vertex v is called the degree of v. Because we only consider 
connected graphs, there are no vertices of degree 0. If a vertex v has degree 1, we call it a boundary 
vertex, otherwise we call it internal. 

It will sometimes be convenient to talk about directed edges of the graph. Each non-directed edge 
produces two directed edges going in the opposite directions. These directed edges are reversals 
of each other. The notation for the reversal of d is d; the operation of reversal is reflexive: d = d. 
Directed edges always come in pairs, in other words, there are no edges that are going in one direction 
only. The set of all directed edges will be denoted by (E. If an edge d emanates from a vertex v, we 
express it by writing v -< d. 

2 Thus for a "typical" eigenvector the notions of "strong" and "weak" nodal domains (see [5]) coincide 
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The number of vertices is denoted by |V| and the number of non-directed edges is \£\. Correspond- 
ingly, the number of directed edges is |(£| = 2\8\. 

Another key definition we will need is the dimension of the cycle space of G. 

Definition 2.2. The dimension i of the cycle space of G is the number of edges that have to be 
removed from £ (leaving V as it is) to turn G into a connected tree. 
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Figure 1. An example of a graph with i = 2. For example, one can cut edges (2, 3) 
and (4, 5) to make it a tree. If is positive on shaded vertices and negative on white 
vertices, the nodal domain count on the graph is 3. On the tree obtained by deleting 
(2, 3) and (4, 5) the nodal count would be 5. 



Remark 2.3. An alternative characterization of i would be the rank of the fundamental group of G. 
There is also an explicit expression for i in terms of the number of edges and number of vertices of 
the graph, 

(2.1) £=\£\-\V\ + l. 
Obviously, i = if and only if G is a tree. 

2.2. Functions on discrete graphs. The functions on G are the functions from the vertex set V 
to the set of reals, xf) : V — > R. We only consider finite graphs, therefore the set of all functions t/j 
can be associated with R' v ', where |V| is the number of vertices of the graph. 

Given a function -0 on G, we define a positive domain on G with respect to if) to be a maximal 
connected subgraph S of G such that if) is positive on the vertices of S. Similarly we define the 
negative domains. Then the nodal domain count z^o ("*/>) is the total number of positive and negative 
domains on G with respect to if), see Fig. [1] for an example. When the choice of the graph is obvious, 
we will drop the subscript G. 

Our interest lies with the nodal domain counts of the eigenvectors of (discrete) Schrodinger oper- 
ators on graphs. We define the Schrodinger operator with the potential q : V —>■ R by 

(2.2) (Ft/,) u = -X> 

The eigenproblem for the operator H is Htf) = Xif). The operator H has |V| eigenvalues, which we 
number in increasing order, 

Ai < A2 < • • • < A|y|. 

This induces a numbering of the eigenvectors: H^ n ' = X n ^ n \ This numbering is well-defined if 
there are no degeneracies in the spectrum, i.e. Xj ^ Xk whenever j ^ k. By fjf(A„) we denote the 
nodal domain count of the n-th eigenvector i/r^ of an operator H . 
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2.3. Functions on metric graphs. A metric graph is a pair (G, {L e }), where L e is the length of 
the edge e 6 £ . The lengths of the two directed edges corresponding to e are also equal to L e . In 
particular, Ld = L-g. 

We would like to consider functions living on the edges of the graph. To do it we identify each 
directed edge d with the interval [0, L d }. This gives us a local variable Xd on the edge which can 
be interpreted geometrically as the distance from the initial vertex. Note that if the edge d is the 
reverse of the edge d then and Ld — Xd refer to the same point. Now one can define a function on 
an edge and, therefore, define a function ij> on the whole graph as a collection of functions {ipd}de<£ 
on all edges of the graph. To ensure that the function is well defined we impose the condition 
ipd{xd) = ipd(Ld — Xd) for all d e <E. The scalar product of two square integrable functions if> and <f> 
is defined as 

rL e 

(2.3) (l/),4>) = ^2 1pe{Xe)M X e)dx e . 

This scalar product defines the space L 2 (G). 

To introduce the main object of our study, the nodal domains, on metric graphs we need to define 
the notion of the metric subgraph of (G, {L e }). 

Definition 2.4. A metric subgraph of (G, {L e }) is a metric graph obtainable from G by (a) cutting 
some of the edges of G and thus introducing new boundary vertices, (b) removing some of the edges 
and (c) removing all vertices of degree 0. 

An example of a metric subgraph is shown on Fig. [2j Now, similarly to the discrete case, we can 
define the nodal count for a real-valued function (p. 




Figure 2. An example of a graph and its metric subgraph. The shaded vertices are 
the new ones which appeared due to cuts. 

A positive (negative) domain with respect to a real- valued function is a maximal connected 
metric subgraph on whose edges and internal vertices (f> is positive (corresp. negative). The total 
number of positive and negative domains will be called the nodal count of 4> and denoted by u(<pi). 

We are interested in the nodal counts of the eigenfunctions of the Laplacian A = — -py. As its 
domain we take the set of continuous function^ that belong to the Sobolev space H 2 (e) on each 
edge e and satisfy the Kirchhoff condition 

(2.4) V -^d(O) = for all v E V. 

^— ' ax 

d>~v 

Note that the sum is taken over all directed edges that originate from the vertex v and the derivative 
(which depends on the direction of the edge) is taken in the outward direction. The Laplacian can 
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also be denned via the quadratic form 

(2.5) Q A [^] = J2 r mx e )\ 2 dx e . 

The domain of this form is the Sobolev space 

For boundary vertices condition (12. 4p reduces to the Neumann condition ip' d {0) = 0. We also 
consider other homogeneous conditions on the vertex v, of the general form 

(2.6) ^d(O) cos a v = ip d (0) sin a v , 

where the Neumann condition corresponds to the choice = 0. The corresponding quadratic form 
will then chang^ to 

(2.7) Q A [<0] = J2f e We(xe)\ 2 dx e + Y, V>»tana v , 

eeS J ° u:deg(u)=l 

where the sum is over the boundary vertices and ip(v) is the value of the function at the vertex v. 

Our results will also apply to Schrodinger operators H = A + q(x) with a potential q(x) which is 
continuous^ on every edge of the graph. 

Schrodinger operator H, defined in the above fashion, has an infinite discrete spectrum with no 
accumulation points. As in the discrete case, we number the eigenvalues in the increasing order. We 
will denote by if)^ the eigenvector corresponding to the eigenvalue A n . 

2.4. Our assumptions and results. Let A n be the n-th eigenvalue of the Schrodinger operator H 
on either discrete or metric graph. Let i/* be the corresponding eigenfunction. We shall make the 
following assumptions. 

Assumption 1. The eigenvalue X n is simple and the corresponding eigenvector if)^ is non-zero on 
each vertex. 

Remark 2.5. The properties described in the Assumption are generic and stable with respect to 
a perturbation. Relevant perturbations include changing the potential {q v } in the discrete case 
and changing lengths {L e } in the metric case. More precisely, in the finite-dimensional space of 
all potentials (corresp. lengths) the set A n on which (A ra , ifi( n >) satisfy the Assumption is open and 
dense unless the graph is a circle (see [11J, where this question is discussed for metric graphs). We 
also mention that on each connected component of the set A n the nodal count of tj}^ 1 ' remains the 
same. Indeed, on discrete graphs the sign of the eigenvector on each vertex must remain unchanged. 
On metric graphs the zeros cannot pass through the vertices. Moreover zeros cannot undergo a 
bifurcation (i.e. appear or disappear) — otherwise at the bifurcation point the eigenfunction and its 
derivative are both zero. By uniqueness theorem for Hi[) = this would mean that ij) is identically 
zero on the whole edge, contradicting the Assumption. 

Now we are ready to state the main theorem which applies to both discrete and metric graphs. 

Theorem 2.6. Let X n and tjj^ 1 ' be the n-th eigenvalue and the corresponding eigenvector of the 
Schrodinger operator H on either discrete or metric graph G. // (A n , i//™- 1 ) satisfy Assumption^ 
then the nodal domain count of is bounded by 

(2.8) n-i < z/(V> (n) ) < n, 

4 if cos a v — — the Dirichlet case — the condition ijj{v) = should instead be introduced directly into the domain 
of Qa 

5 Or has finitely many jumps: the jumps can be thought of as "dummy" vertices of degree 2 
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where £ = \£\ — |V| + 1 is the dimension of the cycle space of G. In particular, when G is a tree, 
utyW) = n. 

While we state the theorem in the most complete form, we will prove only those parts of it that we 
believe to be new. The upper bound on the number of nodal domains is a result with a long history 
going back to Courant P, [7] . The original proof for domains in IR d was adapted to metric graphs by 
Gnutzmann, Weber and Smilansky [T5] . who used the M. d proof from Pleijel [21] who, in turn, cites 
Herrmann [T7] who simplified the original proof of Courant [6] . 

The history of the discrete version of Courant 's upper bound is more complicated. The question 
was considered by Colin de Verdiere [5], Friedman [12], Duval and Reiner [9], and Davies, Gladwell, 
Leydold and Stadler [8]. The latter paper contains a good overview of the history of the result and 
points out various shortcomings in the preceding papers. The point of difficulty was counting the 
nodal domains if an eigenvalue is degenerate (and therefore there is an eigenvector which is zero on 
some vertices). As shown in [8], the upper bound is n + m — 1, where m is the multiplicity of the 
eigenvalue. In our case, Assumption [fl which is essential for the lower bound (see Section [3^1) . also 
simplifies the upper bound. 

The lower bound for the nodal domains on metric trees (i.e. the £ = case) was shown by Al- 
Obeid, Pokornyi and Pryadiev [H [23], [22] and by Schapotschnikow [21] . For completeness, we give a 
sketch of the proof of this case in the Appendix. 

Finally, the results on the lower bound for discrete graphs (both £ = and £ > cases) and for 
metric graphs with £ > are new and will be proved in this paper. 

Note added in proof: It has been brought to the author's attention by J. Leydold that the lower 
bound for discrete trees has been also obtained by Biyikoglu [3] as a corollary of a result of Fiedler 

ma- 



3. Proofs 

We will apply induction on £ to deduce the statement for metric graphs. The proofs for the discrete 
case follow the same ideas but differ in some significant detail. 

First, however, we discuss an important consequence of Remark 12.51 it is sufficient to prove 
statements on nodal counts under the following stronger Assumption. 

Assumption 2. Assumption^ is satisfied for all eigenpairs (A^, xj)^) with k <n. 

Indeed, if only Assumption [1] is satisfied but Assumption [2] is not, we can perturb the problem 
so that (a) the nodal count of the n-th eigenfunction t/r n ) does not change and (b) Assumption [1] 
becomes satisfied for all k < n. Then, anything proved about the nodal domains of ip^ in the 
perturbed problem (which satisfies Assumption [2]) will still be valid for the unperturbed one. 

In our proofs we use the classical ideas of mini-max characterization of the eigenvalues. Let H 
be a self-adjoint operator with domain T>. Assume the spectrum of H is discrete and bounded from 
below. Let Qh[4A = (ij>,H'4 ) ) be the corresponding quadratic form. Then the eigenvalues of H can 
be obtained as 

3.1 X k+ i = max mm — — — , 

h,...,f k ev> vex>, fjW)=o (il>,i/>) 

where the maximum is taken over all linear functionals over T>. 

We will need the following classical theorem (see, e.g., Chapter VI] or [T6J Chapter II]) 

Theorem 3.1 (Rayleigh's Theorem of Constraint). Let H be a self-adjoint operator defined on T>. 
If H is restricted to a subdomain V R = {ip G V : g(if)) = 0}, where g £ V , then the eigenvalues fi n 
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of the restricted operator satisfy 

An < f-n < A n +i, 

where A n are the eigenvalues of the unrestricted operator. 

3.1. Metric graphs (£ > 0). We will derive the lower bound for graphs with cycles by cutting the 
cycles and using the lower bound for trees. 

Proof of Theorem \2.6i for metric graphs (£ > 0). We are given an eigenpair (A n , t/j^). Assume that 
cutting the edges e\, . . . , &t turns the graph G into a tree. We cut each of these edges at a point 
Xj G Cj such that ^ n \xj) ^ 0. We thus obtain a tree with \£ (G)| + £ edges and |V(G)| + 2£ vertices. 
Denote this tree by T. There is a natural mapping from the functions on the graph G to the functions 
on the tree T. In particular, we can think of rj)^ n ' as living on the tree. 

We would like to consider the same eigenproblem Hip = fit/) on the tree now. The vertex conditions 
on the vertices common to T and G will be inherited from the eigenproblem on G. But we need to 
choose the boundary conditions at the 2£ new vertices. Each cut-point Xj gives rise to two vertices, 
which we will denote by Uj + and Define 



U ; 



"?+ ,(n)/ \ ' a 3- i(n)( \ ' J 1 ' 

where the derivatives are taken in the inward direction on the corresponding edges of T. Since ip^ n \ 
as an eigenfunction, was continuously differentiable and i/' (ri '' ) (' u j+) = ^ n \ u j-)> we have 
Now we set the boundary conditions on the new vertices of T to be 

—il)(u j+ ) = a j+ ip(u j+ ), —ip(uj-) = aj-^Uj-), j = 1, . . . , £, 

where the derivatives, as before, are taken inwards. By definition of the coefficients a,j±, the function 
xj}^ 1 ' satisfies the above boundary conditions. It also satisfies the equation Hip = fj,ip and the vertex 
conditions throughout the rest of the tree. Thus, ijy n ' is a ls° an eigenfunction on T and A n is the 
corresponding eigenvalue. If we denote the ordered eigenvalues of T by fik, then A n = fi m for some 
m. It is important to note that m is in general different from n. We will now show that m > n. 

Denote by Q&[ip] the quadratic form corresponding to the eigenvalue problem on G; its domain 
we denote by T~Lq. Similarly we define Qt[i^} and Ti,j. As we mentioned earlier, there is a natural 
embedding of TC& into Ti,j. Moreover, we can say that 

T-Lg = {tp e Hj : if>(u j+ ) = tp(uj-),j = 1,. . .,£}. 

We also note that, formally, 

t 

If tp £ 7i<Q then ip(uj + ) = ip(ujJ) and Oj+ = —aj- result in the cancellation of the sum on the 
right-hand side. This means that on Hq, Qj[tp] = Qg[iJ)]. 

Now we employ the minimax formulation for the eigenvalues A& on G, 

A fc+1 = max min Qg^P] — max m i n Qt[iP], 

\\lp\\=l, 1p±4>i IIV'INL V'-L'/'i 

Comparing it with the corresponding formula for the eigenvalues on T 

fi k+1 = max min Qj[ip], 

IIV»II=1. V»-L<fc 
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we see that the eigenvalues A& correspond to the same minimax problem as fi^ but with £ additional 
constraints = tp(uj-). By Rayleigh's theorem we conclude that Hk < Xk for any k. Therefore, 

if A n = fi m for some n and m, they must satisfy m > n. 

To finish the proof we need to count the number of nodal domains on G and on T with respect 
to V (n) - When we cut an edge of G, we increase the number of nodal domains by at most onej. 
Therefore, 

ur(ijj) < v g (i/j) + t. 

On the other hand, we know that the nodal counting on the tree is exact, and, since ifi^ is the m-th 
eigenvector on T, 

z/ T (i/> (n) ) =m>n. 
Combining the above inequalities we obtain the desired bound 

v G {il> {n) )>n-l. 

To conclude the proof we acknowledge that we implicitly assumed that the tree T satisfies As- 
sumption (TJ more precisely, that the eigenvalue /i r is simple. To justify it, we observe that, if this is 
not the case, a small perturbation in the lengths of the edges will force T to become generic but will 
not affect the properties of the eigenvectors of G. □ 

3.2. Discrete trees (£. = 0). Take an arbitrary vertex of T and designate it as root, denoted r. The 
tree with a root induces partial ordering on the vertices V: we say that v± < v 2 if the unique path 
connecting V\ with r passes through v 2 (see Definition 12. ip . We denote by Vi -< v 2 the situation when 
vi < v 2 and v± ~ v 2 - 

In the above ordering the root is higher than any other vertex. Since T is a tree, for each vertex v, 
other than the root, there is a unique u such that v -< u. Given a non-vanishing i\) we introduce the 
new variables R v = ip u /4>v, where v -< u. Variables R v are sometimes called Riccati variables [20] . 

The eigenvalue condition Hip = \if) can now be written as 

(3.2) - ip u - 22 + qvi>v = tyv, 

w~<v 

and, after dividing by ip v , 

(3.3) R v = q v - TT- 
If v is the root, condition (I3.2p takes the form 

Therefore, if we define 

R r = q r — A — ^ — — , 

then the zeros of R r in terms of A are the eigenvalues of H. Whenever R r (X) = 0, the values of R v , 
«/r, uniquely specify the corresponding eigenvector rj) of H, and vice versa. 

Equation ( 13.31) provides a recursive algorithm for calculating R v , in order of increasing v. Thus 
one gets a closed formula for R v in terms of q u , u < v and A. This is best illustrated by an example. 



'The number of nodal domains might not increase at all if a nodal domain entirely covers a loop of G 
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Figure 3. An example of a tree graph with 5 being the root vertex. 

Example 3.2. For the tree shown in Fig. [3] the eigenvalue condition in terms of Riccati variables 
reads 

(3.4) 

(3.5) 

(3.6) 

(3-7) ■ • ; - •/ 1 • 

(3.8) = q 5 -X- — -^ = R 5 

By substituting lines (13.4p and (13.51) into equation (13.71) . and then lines (I3.6P and (13.71) into equation 
(13.81) . one obtains an eigenvalue condition for H. 

Denote by P v the set of all poles of R v with respect to A and by Z v the set of all zeros of R v ; these 
sets are finite. We define N< to be the number of negative values among R u with u < v ; we similarly 
define N^: 

(3.9) iV< = {u < v : Ru < 0} , = {u < v: R u < 0} 

The above numbers are not defined whenever one of R u has a zero or a pole. The following lemma, 
listing properties of the Riccati variables, their poles and zeros, amounts to the proof of Theorem 12.61 
when G is a tree and q is generic. 

Lemma 3.3. Assume that, for each v, the sets Z w with w -< v are pairwise disjoint for all v. Then 

(1) Pv = \J v ^ v Z u 

(2) For every p G P v , Umx^ p -oR v = —00 and lim\^ p+0 R v = +00. Also, limx^_ 00 R v = +00 
and lim\^ 00 R v = — 00. Outside the poles, R v is continuous and monotonically decreasing as 
a function of A . 

(3) There is exactly one zero of R v strictly between each pair of consecutive points from the set 
{-00} U {00} U P v . 

(4) Between each pair of consecutive points from {—00} U {00} U Z v , the number (where 
defined) remains constant. When a zero of R v is crossed, increases by one. 

(5) Between each pair of consecutive points from { — 00} U {00} U P v , the number N£ (where 
defined) remains constant. When a pole of R v is crossed, N< increases by one. 

(6) When A = A n is an eigenvalue of H , the number of the nodal domains of ip^ is given by 

(3.10) u(X n ) = N< + l. 

Proof. Part [1] follows directly from equation (13. 3p . 

Part [2] follows from (13. 3p by induction over increasing v. 

Part [3] follows from part [2j between each pair of consecutive points from {—00} U {00} U P v , the 
function R v decreases from +00 to —00. 
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Parts H] and [5] are linked together in an induction over increasing v. The induction is initialized by 
for minimal v (i.e. there is no w with w < v). In this case, R v = q v — A, therefore N^- = to the 
left of A = q v and = 1 to the right of A = q v . 

The inductive step starts with part [51 For a vertex v, let both statements be verified for all w, 
w < v . The statement for iV< is obtained immediately from the duality between the zeros and the 
poles (part[Tj). Note that the assumption of the lemma implies that only one of with w -< v can 
increase when A crosses a pole of R v . 

To obtain the statement for consider two consequent poles and two consequent zeros of R v , 
interlacing as follows 

p 1 < z 1 <p 2 < z 2 . 

Then R v is positive for A G (pi, z\) (by part [2]), therefore, on this interval = N<. When z\ is 
crossed, increases by one since R v becomes negative: = N< + 1 = C. On the other hand, 
when p 2 is crossed, and N< become equal again since R v > 0. However, N< has increased by 
one (by the induction hypothesis) and therefore is still equal to C. The above is obviously valid 
even if p% = — oo or/and z 2 = +oo. 

Finally, to show part Owe observe that N< is the number of negative Riccati variables throughout 
the tree. If R v < then the signs of ip v and ip u (where u is the unique vertex satisfying v -< u) are 
different, i.e. the edge (u, v) is a boundary between a positive and a negative domain. Removing 
all boundary edges separates the tree into subtrees corresponding to the positive/negative domains. 
But removing N< edges from a tree breaks it into N< + 1 disconnected components, therefore the 
number of domains on a tree is equal to iVj5 + 1. □ 

Proof of Theorem \2.6[ for discrete trees (£ = 0). The condition of Lemma [3.31 is satisfied due to the 
genericity assumption. Indeed, if there are v, W\ and w 2 such that wi -< v, w 2 -< v and A G Z Wl fl Z W2 
then one can construct an eigenvector with eigenvalue A and with ip v = 0. 

Since the sets Z v are finite, N< must become zero when A — > — oo. Consequently, N< is zero 
between — oo and the first pole of R r . Denote by p n the n-th pole of R r . By part [3] of the lemma, 
the first eigenvalue Ai of H lies in the interval (oo,pi), on which N< is zero. By (13.101) we thus have 
z/(Ai) = 1. Further, A2 lies in the interval {pi,p 2 )- By part [5] of the lemma, N< = 1 on this interval, 
giving u(X 2 ) = 2. Equality for other A n follows similarly. □ 



3.3. Discrete graphs (£ > 0). In this case if is a matrix and the quadratic form is 

|v| 

(3.11) Qo[V>] = QhW\ = H j^k, 

j,k=i 

where 

{-1, 3 ~ k, 
Qj, J = k , 
0, otherwise. 

Proof of Theorem \2.6[ for discrete graphs (£ > 0). We will prove the result by induction. The initial 
inductive step £ = is already proven in Section 13.21 

Assume, without loss of generality, that we can delete the edge (1, 2) of the graph G without 
disconnecting it. We will denote thus obtained graph by T. Note that V(G) = V(T). Let <f> be an 
eigenvector of Hq with eigenvalue A n . We would like to prove that vq{4>) >n — I. 
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Set a = 02/01 and define the potential p on T by 

{Qi - a, j = 1, 
g 2 - 1/a, J = 2, 

With the aid of potential p we define the operator H r in the usual way, see equation (12.21) . It is easy 
to see that, due to our choice of potential p, the vector (j> is an eigenvector of Hr- For example, 

(H r 4>) 1 = - ^ <Pj + (<?i ~ a )0i = -^2<fij ~ 4*2 + gi0i = (H G (f)) 1 = A n 0i, 

where the adjacency is taken with respect to the graph T. 

The eigenvalue corresponding to remains unchanged. However, in the spectrum of Hr, 

this eigenvalue may occupy a position other than the n-th. We denote by m the new position of A n : 

Now consider the quadratic form associated with Hr- Consulting (13. lip we conclude 

(3.12) Qr#] = Qg[V] + 2^2 - oaj>l - a~V 2 - 
Consider first the case a > 0. We write Qr[V'] m the form 

Qr#] = QgM - (« 1/2 Vi - «- 1/2 ^ 2 ) 2 < QgM- 

From here and equation ( 13. lh we immediately conclude that fij < Xj. Therefore, \i m = X n implies 
m > n. 

From the inductive hypothesis we know that vr[<f>] > m—(l — l). But the number of nodal domains 
of <p with respect to F is either the same or one more than the number with respect to G: a > 0, 
therefore 0i and 2 are of the same sign and we may have cut one domain in two by deleting the 
edge (1,2). In particular, z/r[</>] < Ug[4>] + 1. Eliminating [</>], we obtain v&[(f>] + 1 > m — (I — 1), 
which is the sought conclusion. 

In the case a < the quadratic form on T can be written as 

(3.13) QrM = QgM + (/9 1/2 ^i + Z?" 172 ^) 2 , 
where j3 — — a. Consider the subspace 

V R = {^E M |v| : /? 1/2 ^i + P~ 1/2 iJ2 = 0}. 

The restrictions of Hq and Hr to this subspace coincide, as can be seen from (13.131) . Therefore we 
can apply Theorem 13.11 twice, obtaining 

Xj-i < Pj-i < Aj A'i-i ^ Pj-i ^ A*j j 

where pj are the eigenvalues of the restricted operator. In particular, we conclude that jij-i < Xj. 
Since A-spectrum is non-degenerate, < Xj + i, therefore // m = A n implies m > n — 1. 

On the other hand, the number of nodal domains with respect to T is the same as with respect to 
G: since a = 2 /0i < 0, we have cut an edge between two domains. Using the inductive hypothesis 
we conclude that 

u G (<t>) = v T {<t>) >m-(l-l)>n-l-(l-l)=n-l. 

We finish the proof with a remark similar to the final statement of the proof for metric graphs. 
If the new graph F happens not to satisfy Assumption (TJ a small perturbation in q will force F to 
become generic but will not affect the properties of the eigenvectors of G. □ 
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3.4. Low nodal count in a non-generic case. In this section we show that the genericity assump- 
tion (Assumption [Q) is essential for the existence of the lower bound. We shall construct an example 
in which the assumption is violated and the nodal count becomes very low. The construction is 
based on the fact that an eigenfunction of a graph (as opposed to a connected domain in M d ) may 
be identically zero on a large set. 

We consider a metric star graph, which is a tree with N edges all connected to a single vertex. 
For Dirichlet boundary conditions one can show [18] that k 2 is an eigenvalue of the graph if 



To obtain all eigenvalues of the star graph, one needs to add to the solutions of (13.141) the points 
which are "multiple" poles of the left-hand side of (I3.14p . More precisely, if a given A; is a pole for 
m cotangents at the same time, then k 2 is an eigenvalue of multiplicity m — 1. Those eigenvalues 
that are not poles (but zeros) of the left-hand side of (13. 14ft interlace the poles: between each pair of 
consecutive poles (coming from different cotangents) there is exactly one zero. 

Now we choose the lengths Lj to exploit the above features. Let L\ = \,L% = 1/m for some m G N, 
and the remaining lengths be irrational pairwise incommensurate numbers slightly greater than 1. 
By construction, k = mir is a pole for cot(fcLi) and cot(kL2). The corresponding eigenfunction is a 
sine-wave on the edges 1 and 2 and is zero on the other edges. It is easy to see that it has m + 1 
nodal domains. On the other hand, counting the poles of (13.141) . one can deduce that there are 
(m — 1)(N — 1) + 1 eigenvalues preceding (mir) 2 . Thus, we have constructed an eigenfunction which 
is very high in the spectrum but has low number of nodal domains. 

A similar construction is possible for discrete graphs as well. 



The result of the present article came about because of two factors. The first was the request 
by Uzy Smilansky that the author give a talk on the results of [24J at the workshop "Nodal Week 
2006" at Weizmann Institute of Science. The second was the discussion the author had with Rami 
Band on his proof that the nodal count resolves isospectrality of two graphs, one with I = and 
the other with I = 1 (now a part of [2]). Rami showed that the nodal count of the latter graph is 
v(ij)^) = n — 1 or n with equal frequency. His result lead the author to conjecture that for the 
graphs close to trees the nodal count of the n-th eigenstate does not stray far from n. The author is 
also grateful to Uzy Smilansky and Rami Band for patiently listening to the reports on the progress 
made in the proof of the conjecture and carefully checking the draft of the manuscript. 

The author is indebted to Leonid Friedlander for his explanations of the results and techniques of 
[IT] . The author is also grateful to Tsvi Tlusty for pointing out reference [8], to Vsevolod Chernyshev 
for pointing out [2"3"j |2"2"]. to Vladimir Pryadiev for pointing out [I] and to Philipp Schapotschnikow 
for several useful comments. 

Most of the work was done during the author's visit to the Department of Physics of Complex 
Systems, Weizmann Institute of Science, Israel. 



In this section we give an informal overview of the proof of (12.81) on a metric tree (I — 0). For 
detailed and rigorous proofs we refer the reader to [23], [221 124"] . 

Let (A n , if)^) be an eigenpair for a tree T satisfying Assumption [2j Choose an arbitrary boundary 
vertex of the tree T and call it the root r. We can now orient all edges of the tree towards the root 
(well-defined because it is a tree) and will be taking derivatives in this direction. For each non-root 



N 



(3.14) 
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vertex v there is only one adjacent edge that is directed away from it. We call it the outgoing edge 
of the vertex v. The other adjacent edges are correspondingly incoming. An incoming subtree of 
vertex v is defined recursively as the union of an incoming edge (u, v) with all incoming subtrees of 
the vertex u, see Fig. HI 




orientation 
3^ 



Figure 4. An example of a tree with root r. If v is vertex 6 then it has two incoming 
subtrees, one consisting of edges (1,3), (2,3) and (3,6) (highlighted in thicker lines) 
and the other consisting of only one edge (4,6). The outgoing edge of v is the edge 
(6,7). 



If we drop the boundary condition at the root, then for any A < A„ there is a solution <p(X,x) 
which solves the equation H<p = A(/> and satisfies all remaining vertex conditions. This solution is 
unique up to a multiplicative constant. 

The function can be constructed recursively. We fix A and initialize the recursion by solving 
the equation H(f> = \<p on the outgoing edge of each non-root boundary vertex and imposing the 
boundary condition corresponding to this vertex. 

Now let v be a vertex such that the equation is solved on each incoming subtree T v -. We denote 
these solutions (which are defined up to a multiplicative constant) by Cj<pj(x). We would like to 
match these solutions and to extend them to the outgoing edge of v. 

Denoting the solution of the outgoing edge by <p v (x) we write out the matching conditions at the 
vertex v, 

Mv) = C^v) = C 2 <p 2 {v) = . . . 
( j ) ' v (v) = CMv) + C 2( f>' 2 (v) + ... 

Suppose that all of the functions <pj(x) assume non-zero values on the vertex v. Then the condition 
on <fi v (x) takes the form 

\M V ) M v ) 

It is now clear that 4> v , as a solution of H<p = \<fi satisfying this condition, is also defined up to a 
multiplicative constant, C v . The continuity condition now fixes the constants Cj to be C v <fi v (v ) / 4>j(v). 
Thus we obtain the solution on the union of subtrees T v - and the outgoing edge of v. This union is 
in turn an incoming subtree for another vertex (or the root). 

In the case when one of <pj(x) is zero on the vertex v (without loss of generality we take 4>i(v) = 0), 
the condition on <p v takes the form <p v (v ) = 0. The solution <p v is again defined up to a multiplicative 
constant C v . The values of the other constants are now given by C\ = C v (j)' v (v) / 4>'i(v) and Cj = 
when j > 1. Again the solution on the union of subtrees and the outgoing edge of v is obtained 
up to a constant. 

Finally, if more than one of 4>j(x) is zero on the vertex v (without loss of generality, <pi(v ) = <p 2 (v) = 
0), one can take Cj = for all j > 2, find non-zero C\ and C 2 such that Ci^iy) + C2<p' 2 {v) = and 
extend the function by zero on the rest of the tree. This function will satisfy the Kirchhoff condition 
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at v and also all other vertex conditions. Thus it is an eigenfunction and, moreover, it is equal to 
zero at an inner vertex. This contradicts our assumptions. 

We have now constructed a function 0(A, x) which coincides with the eigenfunction of the tree 
whenever it satisfies the boundary condition at the root. To count the nodal domains we need to 
understand the behavior of zeros of <p as we change A. In order to do that we consider the function^ 
R(X,x) = <fi' (A, x) I <fi(X, x) where the derivative is taken with respect to x in the direction towards 
the root. If x is a zero of (p, it becomes a pole of R(X, x). From the definition of -R(A, x) we see that 
R(X, x — 0) = — oo and -R(A, x + 0) = oo. Differentiating R(X, x) with respect to x and using the 
equation —</)" + q(x)(f) = X<p, we see that R(X,x) satisfies 

4-R = q(x) - A - R 2 , 
ax 

a Riccati-type equation. Conditions (12. 6p on the boundary vertices in terms of R(X, x) take the form 
R(X,v) = tan(a„) G M 1 U {oo}- The matching conditions on the internal vertices imply that the 
value of R(X,v) on the outgoing edge is equal to the sum of the values of R(X,v) on the incoming 
edges (in general, R is not continuous on internal vertices). 

Now let A2 > Ai and R(X2,xo) = R(Xi,xo). Then R'(X2,x) < R'(Xi,x) and therefore, on some 
interval (x ,Xo + e), we have i?(A 2 ,x) < R(Xi,x). Moreover, once R(X 2 ,x) < R(Xi,x), we have 
R(X2,y) < R(Xi,y) for all y > x provided both functions do not have poles on [x,y]. This can be 
seen by assuming the contrary and considering the point z G [x, y] where R(X2,z) = R(Xi,z). 

Using these properties one can conclude that for each fixed xo, the value -R(A, xo) is decreasing as a 
function of A between the pairs of consecutive poles. A direct consequence of this is that the poles of 
R move in the "negative" direction as the parameter A is increased. The zeros of 0, therefore, move 
in the direction from the root to the leaves. Since q(x) is continuous, zeros of <fi cannot bifurcate on 
the edges, see Remark 12.51 in Section 12.41 

To see that the zeros of 4> do not split when passing through the vertices, assume the contrary 
and consider the reverse picture: A is decreasing. There are at least two subtrees with zeros of </> 
approaching the same vertex v as A approaches some critical value from above. At this critical value 
we thus have two subtrees on which </> has zero at v. But earlier we concluded that this situation 
contradicts our genericity assumption. 

To summarize, as A is increased, new zeros appear at the root and move towards the leaves of 
the tree. The zeros already in the tree do not disappear or increase in number. Now suppose A& is 
an eigenvalue and thus R(Xk,r) = tan(a r ). As we increase A the value of R(X,r) decreases to —00, 
jumps to +00 (when a new zero enters the tree) and then increases to tan(a r ) again. Thus between 
each pair of eigenvalues exactly one new zero enters the tree. And, on a tree, the number of nodal 
domains is equal to the number of internal zeros plus one. 
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